Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 28 May 2010 (MN mjX style file v2.2) 

Using the Topology of Large Scale Structure to constrain Dark 
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o 

(N 



(N 

O 

u 

43 

Oh 
6 

(N 
> 

m 
\o 
en 

in 

o 
o 



Caroline Zunckel 12 , J. Richard Gott III 1 and Ragnhild Lunnan 3 

1 Astrophysics Department, Princeton University, Peyton Hall, 4 Ivy Lane, NJ, 08544, USA 

2 Astrophysics and Cosmology Research Unit, University of Kwazulu-Natal, Westville, Durban 4000, South Africa 

3 Harvard-Smithsonian Center for Astrophysics, 60 Garden St., Cambridge MA 02138, USA 



28 May 2010 



1 INTRODUCTION 
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ABSTRACT 

The use of standard rulers, such as the scale of the Baryonic Acoustic oscillations (BAO), 
has become one of the more powerful techniques employed in cosmology to probe the entity 
driving the accelerating expansion of the Universe. In this paper, the topology of large scale 
structure (LSS) is used as one such standard ruler to study this mysterious 'dark energy'. By 
following the redshift evolution of the clustering of luminous red galaxies (LRGs) as measured 
by their 3D topology (counting structures in the cosmic web), we can chart the expansion rate 
and extract inf ormation about the equation of state of dark energy. Using the technique first 
introduced in dPark & Kinj d2009l) ). we evaluate the constraints that can be achieved using 3D 
topology measurements from next-generation LSS surveys such as the Baryonic Oscillation 
Spectroscopic Survey (BOSS). In conjunction with the information that will be available from 
the Planck satellite, we find a single topology measurement on 3 different scales is capable 
of constraining a single dark energy parameter to within 5% and 10% when dynamics are 
permitted. This offers an alternative use of the data available from redshift surveys and serves 
as a cross-check for BAO studies. 

cIa(z) — (1 + z)~ 1 r(z), which is related to the expansion rate 
of space via the comoving distance 
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Ustier et alj feOOqY ISpergel et al. ( 20061) : lRiesset al. (2006)). This 
provided convincing evidence for the presence of an unidentified 
entity in the Universe which (in the context of General Relativity) 
must act against the gravitational attraction of ordinary matter. In- 
direct yet also compelling evidence came from the missing energy 
density inferred from the discrepancy between the measurements of 
matter density (from direct measurements and observation of the 
Integrated Sachs-Wolfe effect) and the indications of spatial flat- 
ness from the CMB anisotropy spectrum, as well as the level of the 
initial inhomogeneity measured in the CMB compared with large 
scale structure today. Although the presence of this 'dark energy' is 
now well-established, its nature still evades us and characterizing it 
has become one of the most important topics in cosmology today. 
This is evidenced by the large number of experiments which have 
been proposed and designed with this question in mind. 

One such effort is the Baryon Oscillations Spectroscopic Sur- 
vey (BOSS) which plans to map the spatial distribution of lumi- 
nous red galaxies (LRG) and quasars over 10, 000 sq. deg. of the 
sky. The survey hopes to detect the excess of galaxy clustering at 
100 Mpc/h separations left over from the acoustic oscillations in 
the baryon distribution at the time of last scattering. The change in 
the characteristic scale of this phenonemon from the time of the 
CMB to today is encapsulated by the diameter angular distance 



where 



H( 



h \ n m (i + z) 



fix exp 



1 + w(z>) 
1 + z' 



dz' 



(2) 
where Hq is the Hubble parameter today (note that flatness is as- 
sumed), c is the speed of light, Q m is the current matter density and 
Qx is the current density of dark energy. We focussing on mea- 
suring the equation of state w(z) of the dark energy component, 
which describes the ratio of its pressure to its energy density as a 
function of redshift. Because the scale of the oscillations at the time 
of last scattering is measured precisely from the CMB peak mor- 
phology, the BAO scale becomes a standard ruler. BOSS is fore- 
casted to measure (1a(z) to 1% at various redshifts dSchlegel et al.l 
(2009)), placing constraints on the equation of state of dark energy 
w in Eqn. [2] In this paper, we will use the topology of large scale 
structure as another such standard ruler with which to get a han- 
dle on w(z). Being a tracer of the primordial density perturbations, 
LSS provides a record of the initial conditions and so its topology 
has been used extensively to test our current assumptions of the 
earliest epoch. For example, the measured topology of the Sloan 
Sky Digitial Survey (SDSS) traced by LRGs was shown to be con- 
sistent with expectatio ns from a Universe with Gaussian randon 
phase initial conditions iGott et al. I ( 12009a) ). In this paper, we rec- 
ognize that topology is another measure of clustering or the num- 
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ber of structures per unit volume and thus by mapping how galaxy 
clustering per unit volume evolves with redshift, the expansion his- 
tory of the Universe can be studied. The po tential of thi s dataset in 
dark energy studies was first recognized in JPark & Kim ( 2009)), in 
which the genus statistic used to characterize LSS is adapted to map 
r(z) at vario us epochs, thereby probing w(z). Using the method 
presented in ( IPark & Kiml J2009I) ), we build on this analysis with 
the aim of determining how effective the measurement of topology 
from future LSS surveys such as BOSS will be at elucidating the 
nature of dark energy. We find that the topology is indeed effec- 
tive in placing constraints of the average equation of state (at the 
5% level). The constraints do however weaking when a smoothly 
evolving equation of state is considered, with very little informa- 
tion on a second dark energy parameter delivered. In Section|2] the 
genus statistics are reviewed, followed by a discussion of their ap- 
plication in Cosmology in Section[3] Section|4]gives the details of 
the analysis. In Section[5] the method is applied to the BOSS data as 
well as information about the matter density and Hubble parameter 
from the Planck satellite. 



2 THE GENUS AND RELATED STATISTICS 

Because the perturbations in the primordial density field on scales 
larger than the correlation scale underwent linear growth, the pat- 
tern of matter overdensities today should reflect the distribution of 
these seeds from which they formed - - high density regions such 
as galaxies and clusters of galaxies, are in fact amplifications of 
the primordial density perturbations. This means that the distribu- 
tion of the structure today on large scales gives us a window to 
the co nditions of the U niverse in a much younger state. As pointed 
out in lGott et al.l (1986). this initial state of the perturbations is not 
recorded only in the pattern of the overdense regions; the cosmic 
underdensities also form part of entire structure and it is the cosmic 
sponge on scales large r than the RMS disp lacement of the matter 
that remains preserved ( IPark & Kiml (|2009)). The topology of LSS 
can therefore be used to directly test predictions of our current the- 
ories describing the initial conditions. To measure the topology, the 
number density distribution is smoothed with a Gaussian smooth- 
ing ball of radius R g . We then find the iso-density contours of the 
smoothed distribution which divides the space into two, where the 
fraction volume on the high density side given by 

f(v) = -±= r e~ x2/2 dx- (3) 



Here v is merely a label for the contours. For example, v — 
denotes the iso-contour of the median density which encloses / = 
50% of the volume, while v — 1.5 labels the / = 7% contour. One 
can then define the genus as a function of u, given by the difference 
between the number of donut-like holes and isolated regions; 



G(u) — No. of holes — No. of isolated regions. 



(4) 



For a random Gaussian phase density field, the genus per unit 
volume, q (v) = G { v)/V, is predict ed from theory to be (c.f. 
iHamiltonetall J1986h ; lGott et al.l Jl987tV ) 



g {v) = A{\~v 2 )e^l 2 
where the amplitude is given by 

«fc 2 >/3) 3/2 
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selected smoothing scale and is independent of the amplitude of 
P(k). Recent studies of the topology of LSS traced by LRGs 
in the SDSS survey iGott et al. I y009bj) at smoothing lengths of 
21 h _1 Mpc and 34 h _1 Mpc show the amplitude of the genus curve, 
A, is modeled very well by N-body simulations and closely fol- 
lows that of the initial conditions. Note that the topology of a non- 
Gaussian density field at any scale within the linear regime is also 
preserved in comoving space. 

The currently favoured model for the beginnings of structure 
formation, called Inflation, assumes that the primordial fluctuations 
to be Gaussian random ph ase, which has bee n shown to lead to a 
sponge-like topology I Gott et al.l dl987l . ll986T) ). The medium den- 
sity contours measured for various LSS samples over the past few 
years have been foun d to be consistent w i th expectations of t his 



type of topology ( see IGott et al. I J2009b1); IPark. eTall d2005al ltah: 
IGott et al. f ( f2009al [i989); Hami lton et all J1986I) ). In addition, the 



Here (k 2 ) is the average value of k 2 in the smoothed power spec- 
trum. A is sensitive to the slope of the power spectrum near the 



theoretical prediction for Gaussian random phase initial conditions 
in Eqn.[5]was shown to provide a suitable fit to the observed genus 
curves for these data sets. In this paper, we focus of whether the 
same data can be used to study the late-time behavior of the Uni- 
verse. 



3 USING TOPOLOGY TO CONSTRAIN COSMOLOGY 

The comoving distance r(z) in Eqn. [T]tells us how scales change as 
a function of redshift due to the expansion of space. This means that 
comparison of characteristic scales in cosmology, such as features 
in the power spectrum or the correlation function, from one epoch 
to the another, provides a handle on r(z). Because the shape of 
the primordial power spectrum P(k) on large scales (where we are 
still in the linear regime) should be conserved in redshift space (the 
transfer function does not depend on z), we can determine P(k) 
for a given set of cosmological parameters from two measurements 
of P(k, z) and expect to find it unchanged. If the power spectrum 
does change, the r(z) relation in[T|and hence the chosen cosmology 
must be wrong. 

As described in Section[2] the genus statistic is related to the 
entire shape of the power spectrum and essentially measures its 
slope near the smoothing scale. The statistic g(v) is a measure of 
the number of structures after smoothing on a certain scale, per 
unit volume, out to a given distance. The measure therefore re- 
lies on a choice of r(z) relation (and hence cosmology) to com- 
pute the enclosed volume. Assuming the incorrect cosmology and 
hence expansion rate will lead to an incorrect estimate of the vol- 
ume in which one is counting structures. In so doing, a different 
amplitude of the genus curve from the true value will be measured. 
Furthermore, the chosen smoothing scale which determines the size 
of the structures will be incorrect. For example, in the case where 
one overestimates the expansion rate, what is assumed to be a unit 
volume will enclose less structure than at z = 0. As a result, the 
measured amplitude of structure on the smoothing scale will be 
lower. At the same time, because the smoothing scale per unit vol- 
ume is kept fixed, smaller structures are effectively smoothed over 
so these two effects partially cancel. Fortunately, there is a net ef- 
fect because the density perturbation field is not scale free and the 
power spectrum has a shape with a varying slope over the region 
of interest. Fig. Q] shows the genus curves for two different cosmo- 
logical models. In this paper, we will exploit the dependence of the 
genus curve on r(z) in order to constrain the equation of state of 
dark energy, as explained in the next section. 

We consider a flat Universe filled with cold dark matter, 
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baryons and dark energy. The cosmological parameters of interest 
are the matter density, Q, m (baryonic and dark), the Hubble param- 
eter Ho and the equation of state of the dark energy component, 
which has been shown to provide a suitable phe nomenological de- 
scription of dark energy JFrieman et alj ( 1200810 . To start, we as- 
sume the simplest parameterization, namely a constant equation of 
state w. The comoving distance along the line of sight, rn (in red- 
shift space), and in the traverse direction, r±, of a feature sitting at 
a redshift z, are related to the redshift range Az covered and the 
angle subtended AO respectively, by 



Az 
Z HJ7) 



c(l + z)d A (z)AQ. 



(7) 



Suppose a length along the line of sight spans a redshift range 
of Az. The comoving distance of the length is given by 
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r \\( z ) 
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If we then perturb the cosmology, the same redshift range Az will 
correspond to a different comoving distance given by 



A,, 



cAz 
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Thus the different comoving distances are related by 

, _ H(z,Q. m ,w,H Q ) 
x ~ x H{z,W m ,w>,H' o y 

The comoving distance subtended by an angle AG is given by 

X z = r-xO) = (1 + z)d A (z, i1 m ,w, Ho) A9 (11) 

If we again perturb the cosmology, the same angle AO will subtend 
a new comoving distance of 



X' z — (1 + z)d A (z,Q' m 
So we can write 



>j',K)AQ 
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(1 + z) d A (z, fim,W, Ho) 
d A (z,Q' m ,w',Ho) 
d A (z,Q m ,w,H ) ' 



(13) 



Similarity, the comoving distances subtended by the angle 
A6 in the y direction in two different cosmologies are related by 



, _ d A (z,n' m ,w',Ho) 
y y d A (z,Q m ,w,H ) 

We can then write the volume of the corresponding ellipsoid 
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so we can relate the comoving volumes in the two different cos- 
mologies via 
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We denote the factor relating the volumes in the different cosmolo- 



gies as 



„ _ d A (z,Q m ,w,Ho) 2 
eff ~ H{z,n m ,w,H ) ' 



(18) 



Re-arranging the above gives an equation which gives the volume 
in the new cosmology as a function of the original volume; 

Vj f . 

(19) 



V' = V 



V eff 
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Suppose we have a galaxy redshift catalogue with which we 
wish to constrain the underlying cosmology. After smoothing the 
distribution with a Gaussian smoothing sphere of dimensions A = 
\x = A y = \ z and volume V = A 3 , we measure an amplitude of 
G — g(X, z)V — g(\, z)X 3 . Now suppose that the assumed cos- 
mology for the smoothing process is in fact incorrect. This means 
that the r(z) relation is wrong. In the true cosmology, the smooth- 
ing sphere is in reality an ellipsoid with volume V' — X' x X' y X' z . 
The same amplitude G = g'(A', z)V is measured by smoothing 
a different cosmology with an ellipsoid or a sphere of radius A', 
which is calculated using 
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(20) 



So we can relate the smoothing scales within the two different 
cosmologies using 



A' = A 



K'// 



V, 
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(21) 



Equating gV = g'V' gives us 
#'(A',2) 



n \ V 



= ff(A,z) 



Vc 



eff 



V iff 



(22) 



Eqn.[22]is very useful because it allows us to calculate the ex- 
pected genus curve for any cosmological model provided we know 
the theoretical g(X,z) at a given scale A for one set of param- 
eter values. This is readily computed using the smoothed matter 
power spectrum g enerated using a Boltzmann code such as CAMB 
JLewis et alj ( 120001) ). The validity of Eqn. |22]depends on how well 
a galaxy distribution smoothed with a sphere of geometric mean 
R g represents that smoothed using an ellipsoid of median radius 
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4 ANALYSIS 

The aim of the analysis is to perform dark energy forecasts for a set 
of measurements of the genus at various smoothing scales, based on 
the 3D distribution of LRGs as it will be measured by BOSS. Using 
a 2048 3 particle cold dark matter N-body simulation. iKim & Parkl 
(2006)'s technique which identifies LRG galaxies by selecting the 
most massive bound halos, corr ectly reproduces the 3D topology 
of the LRG galaxies in the SDSS lGott etaTI J2009bh . 
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We apply the same technique to the Horizon Run Simulation, 
which is the larg est N-body simu lation to date with a volume of 
6.543 (Gpc/h) 3 feimetal] il2009V). Provided we remain in the lin- 
ear regime, the genus curve for the initial conditions remains essen- 
tially unchanged and should provide a suitable representation of the 
genus that would be measured using LRGs as a probe of LSS to- 
day. Hence we use the initial density field from the Horizon Run 
Simulation to construct a hypothetical data set and subdivide the 
simulation into 8 cubes. Placing an observer at the centre of each 
cube, a sphere of comoving radius 1570 Mpc/h (corresponding to 
a redshift of z = 0.6) is carved around the observer and quartered 
to produce 4 BOSS mock catalogues of n steradians and in total, 
32 mock galaxy catalogues with which to test the statistics. The 
genus per unit smoothing volume is then measured for each cata- 
logue at 3 different smoothing scales; 15, 21, 34 Mpc h _1 . These 
scales were selected because they are sufficiently disparate that dif- 
ferent types of structures are being smoothed and can, for statistical 
purposes, be regarded as independent. This gives three independent 
data points gdata(z, As) where i — 1 — N z per redshift bin. The 
noise is taken to be the standard deviation of the measured g values 
of the 32 simulated maps. 

First, we compare the genus per unit smoothing volume mea- 
sured in our mock surveys to the theoretical values from Fig.Q] As 
it should, the shape of the genus curves follows that of Gaussian 
random fluctuations (Eqn. O, but we systematically find slightly 
lower amplitudes than predicted from theory. In particular, the sup- 
pression is systematically larger with larger smoothing length, and 
in redshift bins with more surface pixels. We believe this is due 
to two related effects: pixels on the edges are not counted by the 
program measuring the genus, and structures near the edges are 
supressed by the smoothing and "cut off", both lowering the value 
of the genus measured. To account for this, we introduce correction 
factors for each redshift bin, so that for the original density field, the 
mean values for the given bin correspond to the theoretical value. 

Next, to test the validity of Eqn. [22] we remap the density field 
in each survey, according to two different values of constant w. 
—0.9 and —0.8. For each cosmology, we construct a new density 
field from the original, based on the distances the observers would 
infer according to those cosmologies. We then again compute the 
genus per unit volume for the same three smoothing lengths and 
redshift bins, in our 32 mock surveys. Each bin is corrected for edge 
effects by the factors calculated above, and then compared to the 
expected number based on the geometric mean volume (Eqn. 1221). 
We find that for both cosmologies, in all of the redshifts bins, for 
all smoothing lengths, the differences between the observed values 
in the simulated cosmologies (with w = —0.8 and w = —0.9) 
and theoretical values using Eqn. [22] are well within one stan- 
dard deviation of the mean. We conclude that as long as edge ef- 
fects are taken into account, the theoretical value based on effective 
smoothing length (Eqn. 122b agrees excellently with the observa- 
tions. Therefore, we go ahead and use Eqn. [22] to compute genus 
curve amplitudes in different cosmologies for the likelihood analy- 
sis. 

We now wish to evaluate the likelihood of the parameters 
p — (Ho,n m ,w) in light of the data. This involves computing 
the expected values of the genus per unit volume at the point pi in 
parameter space using Eqn. [22] for a given smoothing length and 
comparing them with the data using 
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Figure 1. Plot of the amplitude of the genus curve gFti as a function of the 
smoothing length R g for the best-fit parameters values for WMAP 3-year 
(n m = 0.24, H = 72) and WMAP 5-year data (Cl m = 0.26, H = 74) 
assuming a flat ACDM cosmology. 
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and Gij is the uncertainty associated with the data smoothed on 
the i th scale in the j" 1 redshift bin. This is computed from variance 
measured from the 32 mock surveys. The likelihood of the trial 
point is then £ (pi\gdata) ■ 

In the above, g re j (A, z) are the expected genus values for our 
reference cosmology, calculated from the linear matter power spec- 
trum generated using CAMB and smoothed on the scale 
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The value of g re f (A', z) at this new smoothing scale A' is found us- 
ing interpolation. We use the WMAP 3-year best-fit cosmological 
parameters to generate g re f (A, z) shown in Fig. [Tj This procedure 
is repeated at each point in parameter space, which is efficiently 
sampled using a Monte Carlo Markov chain (MCMC) algorithm 
until convergence. The 68% and 95% confidence intervals for the 
parameters (Ho, fix , w) are then computed from the chains. 



5 RESULTS 

Assuming a completed BOSS survey, we constructed a hypotheti- 
cal data set consisting of three genus measurements, each at one of 
the three selected smoothing scales. We pick a dark energy model 
with w 7^ —1 to test the ability of the data to distinguish such 
a model from cosmological constant. We simulate a flat Universe 
with a dark energy component described by w — —0.9. 

Fig. [2] shows the 68% and 95% confidence regions in pa- 
rameter space when (a) Ho is fixed while Q, m and w are allowed 
to vary and (b) Q m is fixed while Ho and w are allowed to vary. 
The elongated contours exhibit strong deneracy between the matter 
density and Hubble parameter and the equation of state. For a more 
negative choice of w than the true value, the dark energy compo- 
nent is less important in the past while the matter density today is 
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Figure 2. 95% confidence regions (a) for a parameter space with (O m , w) 
where Hq has been fixed and (b) for a parameter space (w,Hq) where 
Sim has been fixed. The true values of the parameters used to produce the 
simulated topology data are marked by a yellow cross. 



held fixed. In order to match the same genus measurement, a higher 
value of the Hubble constant today is needed in order to restore 
H(z) and to yield the correct measurement of the sampled volume. 
For a more positive choice of w and fixed Ho, H(z) is higher, indi- 
cating that space is expanding more rapidly. The implication is that 
the estimated volume out to a given z will be smaller. To counteract 
this, the matter density can be reduced thereby increasing the con- 
tribution from the dark energy component to restore the expansion 
rate to a lower value. 

In addition to the simulated data from BOSS, we choose 
to include information from the Cosmic Miscrowave Background 
(CMB). The presence of dark energy impacts the CMB primarily 
through the distance to the surface of last scatter, with the expansion 
rate dictating the scales on which the acoustic features appear in the 
CMB spectrum today. The size of the largest full acoustic compres- 
sion which manifests as the first peak in the CMB angular power 
spectrum, is determined solely by the distance that these waves 
could have traveled in the time before recombination, namely the 



sound horizon r s at last scattering; 

rt cmb 

'S \Zcmb, i£b ; i &r J — / Cs&T 

JO 

where c s is the sound speed defined by 



3 1 



3^6 

4^ 



-1/2' 



(26) 



(27) 



Since the baryon to photon ratio can be precisely measured from 
the acoustic peak morphology in the CMB, we can predict c s and 
thus the scale of the first peak. The angular scale on which this 
feature appears today, 6 S depends on the expansion history and the 
distance to the surface of last scattering, d s i 3 ; 
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where 
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(28) 
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Assuming parameter values near the A concordance model, we use 
the fitting formula 

,„ ,o, -0.252 /_. , 2 \ -0.083 

r. = 144.4Mpc ( %#) (%£) . (30) 
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0.14 



This assumes that energy density contribution from the dark energy 
component during this epoch is sufficiently small such that it does 
not drastically affect the size of the acoustic features at the surface 
of last scattering. The information regarding the expansion rate and 
hence dark energy available in the CMB data comes primarily from 
the distance to the last scattering surface d s i a . 

The current constraints from the WMAP 7-year data are 9 S = 
0.597 ± 0.0016 with the matter and baryon densities estimated to 
be Q m h 2 = 0.1098±0.0058and n b h 2 = 0.0250±6.3 x 10~ 4 re- 
spectively (Komatsu et al. (2010)). We also include priors based on 
forecasts for the Planck CMB experiment which should be avail- 
able at the time that the BOSS survey is completed. The higher 
resolution and sensitivity of the Planck Satellite lead to tighter fore- 
casted constraints of 6 S = 0.597±3.1 x 10" 4 , Q. m h 2 = 0.1098± 
1.4 x 10" 3 and Q b h 2 = 0.0250 ± 1.6 x 10~ 4 JColombo et ail 
(2009)). As a result of the degeneracy between w and Ho, permit- 
ting non-A models and allowing w to vary significant degrades the 
constraints on Ho. We thus include the current measure ment of the 
Hubble constant from Ho = 74 ± 3.6km s _1 Mpc _1 (Ri ess et al.l 

Including information from the CMB is highly advantageous 
as the direction of the degeneracy between fl m and w is perpendic- 
ular to that in the topology data. As w becomes less negative, its 
contribution to the overall density at earlier epoch increases. As a 
result, the expansion rate is reduced, with the features on the last 
scattering surface subtending smaller angles, manifesting as a shift 
in the acoustic peaks to lower t. This effect on the CMB anisotropy 
spectrum can however be countered by an increase in the matter 
density which delays recombination and leads to a larger sound 
horizon r 3 . Fig.|4]shows the allowed regions in parameter space for 
the simulated BOSS topology data in conjuction with the WMAP 
7-year (left) and the Planck data (right). The change in the direc- 
tion of the w,£l m ellipse points to degeneracy breaking. There 
is neglible improvement on the dark energy EOS constraint from 



-0.92 4 



to w 



-0.92±° Q °l when the CMB data set is 



updated from WMAP-7 year to the Planck data, despite the clear 
reduction in the confidence interval evident in Fig. [4] This indicates 
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the the topology data is primarily responsible for the constraints on 
w. 

A useful task would be to determining the constraining power 
of this dataset as a function of redshift. This can be studied by 
performing principal component analysis, which essentially iden- 
tifies the directions in the data in which the variation is maximal. 
These directions are captured by z-dependent eigenvectors, with the 
eigenvalue of mode representing how well it is constrained (ie. the 
strength of the information it represents). This will reveal where in 
redshift space this particular dataset is most sensitive to w. 

We start by dividing the redshift region in 10 equally sized 
bins centered at redshifts Zi. We parameterize the EOS as a piece- 
wise constant function, with a constant Wi in each bin where 
Zi — Az/2 < Zi < Zi + Az/2. We then sample the posterior 
probability distribution an MCMC algorithm, marginalizing over 
the other cosmological parameters. Using the chain of samples, we 
construct the covariance matrix as follows; 



C i 



fc=i 



N -1 



(31) 



where X% is the value of the parameter at the i l point in the MCMC 
chain and Xi is the sample mean. We find the eigenvectors by di- 
agonalizing the covariance matrix and decomposing it as follows; 



C = W 1 AW 



(32) 



where the diagonal matrix A contains the eigenvalues correspond- 
ing to the eigenv ectors &j in the rows of matrix W. We follow 
dTang et al.1 (2008)) and track the redshift sensitivity of each mode 
by plotting the quantity 4>i( z )'> 

(33) 



4>i{z) = N | VXiei(z) | 

where the amplitude is proportional to the accuracy with which the 
eigenvector e;(z) can be measured. N is included to lessen the 
dependence on the number of bins. Fig.[3]plots (j>i(z) as a function 
of z. The behaviour of the first and best-constrained mode (red) 
indicates that the constraining power of the survey is a decreasing 
function of z, with most of the information about w delivered at 
low redshift. All other eigenvectors are relatively noisy with no z- 
dependent shape distinct shape. We infer that the redshift sensitivity 
of this data set is likely to behave as <f>\ (z) (shown in red). 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 



Figure 3. Figure showing a plot of tjn (z) as a function of reshift. The first 
eigenvector </>i (z) is shown in red, the second <j>2 (z) in blue and with the 
remaining 8 shown in black. This illustrates the redshift sensitivity of this 
dataset. 





-0.7 -1.3 



(a) 



(h) 



(c) 



5.1 Redshift bins 

Up to this point we have smoothed the entire observed volume on 
three different scales, constructing a dataset comprising of three 
data points. In this Section we consider the effect of dividing the 
observed space into redshift bins of equal volume, thereby increas- 
ing the number of independent measurements three fold. The char- 
acteristics of the redshift shells in the case of one, two and three 
bins are summarized in table l5.ll 

The results of the analysis are summarized in the first section 
of table |5.2| An increase in the number of redshift shells is shown 
to make little difference to the constraints when a constant equation 
of state is considered. The increase in the number of data points is 
clearly countered by the reduction in signal-to-noise in each red- 
shift shell over a single bin. In light of the finding that the sensi- 
tivity of the topology data to the dark energy model is a decreasing 
function of redshift (shown in Fig. [3), it is likely that most of the 
information available about w is concentrated in the first redshift 
bin, with the high-z bins adding little to the analysis. 



(d) 



Figure 4. 95 % confidence regions for (w , Ho ) for the topology data includ- 
ing (a) WMAP 7 year data and (b) Planck data. Plots of (w, O m ) for the 
same combinations of data sets are shown in (c) and (d). The true parameter 
values are shown by the yellow crosses. 



Bins 


Redshift range 


Distance (Mpc/h) 


Volume (Gpc/h) 3 


1 


z < 0.6 


r < 1570 


4.036 


2 


z < 0.46 
0.46 < z < 0.6 


r < 1246 
1246 < r < 1570 


2.016 
1.997 


3 


z < 0.396 

0.396 < z < 0.513 

0.513 < z < 0.6 


r < 1089 
1089 < r < 1371 
1371 < r < 1570 


1.342 
1.326 

1.321 



Table 1. Table summarizing the partitioning of redshift space. 



Using the Topology of Large Scale Structure to constrain Dark Energy 



Number of z bins 




Constant EOS: 


w = -0.9 




1 
2 
3 


-0 Q2+ 005 

u ' a -0.05 

-0 Q2+ 005 
u - yz -0.05 

-0 Q2+ 004 
u - a -0.05 


- 


Dynamical: 


w = -1 


w a = 0.5 


1 

2 
3 


-i.o4±°:i 

-0.93+0-1 

-0.98«-i 


0.45±°;| 

2S +0 ' 5 
041 +0.7 



Table 2. Table summarizing the constraints on the dark energy EOS param- 
eters for different binning schemes using the topology data and the Planck 
data. We assume the EOS to be constant from the last redshift bin at z = 0.6 
to the last scattering surface at z = 1089. 



Figure 5. Plot of the 68% (dark blue) and 95% (light blue) confidence 
regions for (wQ,w a ). The simulated data was divided into 3 equi- volume 
redshift bins summarized in table l5.ll We include the Planck priors as well 
as the current constraint on the Hubble parameter. 



5.2 Constraining dynamical dark energy 

Up to this point, we have selected the simplest description of the 
equation of state, a constant w. Although it has been shown to be a 
good approximation to a quintessense component obeying tracker 
solutions lEf stathioul i 1 9991) ) while accommodating the vacuum en- 
ergy (too = —1) as one of its solutions, it cannot however, be 
used to characterize the effect of scalar field models , in general, or 
modifi ed gravity models on the observable Universe JFrieman et al.l 
(2008)). In this Section, we wish to determine the effectiveness of 
the topology data in constraining a second dark energy parame- 
ter. A popular two-parameter ansatz for w(z ) was introduced in 
dChevallier & Polarskil J200ll) ; iLinded J2003h ) with the following 
form: 



){z) - w +w a (l - a). 



(34) 



where a = 1/(1 + z). This particular function is a favourite in the 
cosmology community because it solves the divergence problem at 
high redshift accompanying other parameterizations, bu t at the cost 
of a m ore rigid assumption of the behavior a priori (Rie ss et al.l 
(2004)). It is also limited i n how well it can co pe with a rapidly 
evolving equation of state JLiddle et al.J |2006)). We simulate a 
set of genus values for a dark energy model with vjq = — 1 and 
W a — 0.5 in three different binning schemes and add Gaussian 
noise. Table |5^2l summarizes the results of the likelihood analysis 
when the simulated topology data in conjunction with the Planck 
priors is used to constrain p = (fl m , Ho,wo, w a )- The constraints 
on wo weaken as expected to roughly 10%. We find that the topol- 
ogy data is relatively uninformative with regards to a second dark 
energy parameter, regardless of the number of redshift bins used. 
In fact, given that data is too noisy to detect any variation in w 
with redshift, partitioning the data into redshift bins serves only to 
weaken the constraints on the average w. Figure[5]shows the 68% 
and 95% confidence regions in (wo,w a ) when the data is binned 
into three redshift shells of equal volume and serves to illustrate 
the increase in the confidence region accompanying the increased 
freedom in the dark energy model. Table |5.2l also indicates that the 
use of a single redshift bin is prefered overall. 



6 DISCUSSION 

The question of the nature of dark energy is one of profound impor- 
tance in cosmology today. The currently favoured model proposes 
that it is energy density associated with the vacuum, with a constant 
EOS of w = —1, identified as the mathematical equivalent of cos- 
mological constant A in the Einstein field equation. In the context 
of our current theories of structure formation, A appears to be most 
successful in reproducing a wide range of present-day observations. 
However, it faces serious theoretical opposition and several recent 
works have highlighted the biases towards w — — 1 that could po- 
tentially arise from t he inclusion of a pri o r information when fit- 
ting the current data teassett et aT] i2004h ; lLmde3 J2004h ). Future 
experiments offer the exciting prospect of determining whether the 
cosmological constant is indeed the correct model. BOSS is fore- 
casted to place constraints on w at the level of a few percent using 
clustering of LRGs on the scale of 100 Mpc. IPark & Kiml b009n 
proposed a way to use the nature of LRG clustering on a wider 
range of scales as measured by the 3D topology to probe the expan- 
sion rate. Because genus statistics relate to the overall shape of the 
power spectrum as opposed to a single scale, they are robust against 
fortuitous noise measurements. Furthermore, the toplogy measure- 
ments offer an alternative use of the data from redshift surveys and 
may be used as a cross check of the conclusions drawn from other 
techniques, such as the BAO method which uses the oscillations in 
the matter power spectrum to extract information. 

In this paper, we have evaluated the constraints on w(z) that 
are achievab le from the BOSS survey using the method presented 
in JPark & K im (2009)). We found that the BOSS survey using the 
topology method alone is capable of placing constraints of 5% on 
a constant equation of state. Allowing for the possibility of time 
variation in w(z) degrades the constraints on wo to 12%, while 
providing weak evidence for a second dark energy parameter. Al- 
though the topology measurement of BOSS may not be capable of 
testing for dynamical behaviour in w(z), it has been shown to pro- 
vide a robust measurement of the average equation of state, which 
may be sufficient to rule out non-A models. 
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